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A dilute atomic-molecular fermion in a trap acquires a 
rather complex layered spatial structure. At very low temper- 
atures, typically at the center of the trap, a superfluid bosonic 
core is formed out of dimers, weakly bound pairs of fermionic 
atoms. Surrounding this core is a mantle, formed by a mixture 
of superfluid bosonic dimers and superfluid fermionic atoms. 
The outermost layer consists of only fermionic atoms, forming 
mostly a normal Fermi gas. 

PACS numbers: 03.75. Ss, 03.75.Hh 

Jin's group [1] has demonstrated that weakly bound 
dimers of fermionic atoms are easily created in a cloud of 
cold Fermi atoms, by ramping up and down the magnetic 
field near a Feshbach resonance. Several other groups 
have been able to quickly confirm these results, by using 
similar experimental techniques [2-4] . In all these exper- 
iments a relatively cold Fermi gas of atoms, interacting 
with a negative scattering length, was driven past the 
Feshbach resonance into a mixture of fermionic atoms 
and bosonic dimers, all of them interacting with positive 
scattering lengths [5] . The ratio of atoms to molecules in 
these experiments could be varied essentially at will, by 
tuning the ramping rate of the ambient magnetic field. It 
was also explicitly demonstrated that this mixed system 
of fermionic atoms and bosonic dimers has a rather long 
lifetime [2,6] and that the rate of various recombination 
processes is rather low, an effect qualitatively suggested 
in Ref. [5] . In a series of related experiments it was shown 
that the weakly bound dimers form a Bose-Einstein con- 
densate [7]. 

The natural question arises: "What is the nature of 
this new quantum object created in these experiments?" 
As it was shown in Ref. [5] , in the regime where the scat- 
tering length is large, but at the same time the various 
average inter-particle separations are even larger, it is 
possible to perform a rather straightforward theoretical 
analysis of its properties. If the system could be regarded 
as locally homogeneous, then by changing the scattering 
length from negative to positive values part of the sys- 
tem undergoes a BCS —* BEC crossover [8-11]. At zero 
temperature, when the scattering length is negative, one 
would expect that the atomic Fermi system would be in 
a BCS state at sufficiently low temperatures. By chang- 
ing adiabatically the scattering length to positive values, 
the size of the Cooper pairs would gradually decrease, 
until dimers, with sizes smaller than the average inter- 



particle separation, are eventually formed and the entire 
system turns into a pure BEC. At the finite temperatures 
achieved so far in atomic Fermi clouds, typically only a 
fraction of the atoms are converted into dimers and one 
ends up with a mixed system of fermions and bosons. 

In all these experiments [1-4] the estimated final tem- 
peratures were well below the critical BEC temperature 
of the Bose component. One could have naively con- 
cluded that a novel mixture of a Bose superfluid mixed 
with a normal Fermi gas was thus realized. One reason 
to believe this would be the fact that the atomic scatter- 
ing length is positive, which corresponds to an effective 
repulsion at these energies. Therefore the Fermi sub- 
system could be expected to be a normal Fermi system. 
The other reason would be the fact that the dimers were 
expected to repel each other [5,11-13] and thus they can 
form a BEC. This aspect was subsequently confirmed ex- 
perimentally [7]. A third reason would be the fact that 
atoms and dimers repel under these conditions too [5,13]. 

Nevertheless, a somewhat subtle quantum phe- 
nomenon could lead to rather dramatic changes of the 
properties of the fermionic component. Namely, even 
though two fermions effectively repel each other (since 
the scattering length is positive), there is a boson in- 
duced interaction, which is attractive and relatively 
long ranged. A fermion traveling in a Bose superfluid 
(fermions always move) will excite a Bogoliubov sound 
mode of the Bose subsystem, which in its turn will prop- 
agate to rather large distances, where it will "shake" an- 
other fermion. This type of interaction between fermions 
is reminiscent of the phonon mediated electron-electron 
interaction in metals, which is responsible for the super- 
conductivity. Technically, this phonon exchange between 
two fermions is a second order perturbation theory con- 
tribution to the ground state energy of the system, and 
as such is always attractive. As a result, if the tem- 
peratures are low enough and the induced attraction 
strong enough, a homogeneous mixture of such bosons 
and fermions could turn into a mixtures of a Bose super- 
fluid in interaction with a Fermi superfluid. If the system 
is dilute and if all the scattering lengths are known [5] and 
the energy of such a system is easily calculated. Here, I 
shall primarily be concerned with the spatial profile and 
spatial structure of such a mixed atom-dimer system at 
low temperatures. Since the induced pairing correlations 
in the Fermi subsystem are in the weak coupling regime 
[5], the presence or absence of the pairing correlations 
plays essentially no role in determining the spatial distri- 
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bution of such an atomic cloud. 

The properties of these novel Bose-Fermi mixtures of 
quantum fluids could be affected by the shape of the 
trapping potential in a non-trivial manner. As it will be 
shown here, the trapping potential has a rather strong in- 
fluence on the detailed spatial matter density proflle and 
spatial composition of these new objects. Even though 
one can create either a pure Bose system or a pure Fermi 
system, however, such objects are likely an exception 
rather than the rule. Most of the time one would cre- 
ate an object, which at the center will be a pure Bose 
superfluid, followed by a layer of mixed Bose and Fermi 
superfluids, with an outer layer being a normal Fermi 
fluid. At the interface between the normal and the super- 
fluid Fermi systems, another pure quantum phenomenon, 
known as Andreev reflection, leads to the formation of 
superfluid correlations in the outer normal Fermi fluid. 
These pairing correlations will be present in a layer of 
width comparable with the so called coherence length 
[14]. When a fermion, impinging on a normal-superfluid 
interface from the normal part, undergoes an Andreev 
reflection, it is reflected essentially backward as a hole. 
This process leads to the appearance in a relatively nar- 
row region of the normal Fermi fluid of pairing correla- 
tions, even though the atom- atom interaction is repul- 
sive. Again, the presence of such pairing correlations 
will have an exceedingly small influence on the spatial 
distribution of the atoms [15]. 

In Ref. [5] it was shown that the energy density of a 
mixed fermion-dimer system in a trap is given by the 
following expression: 
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where V is the trap potential, Ub is the dimer number 
density, n/ = kp/Sn"^ is the total number density of 
both fermion species, assumed here to be present in equal 
amounts, hp is the Fermi wave vector and where 
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and a is the atom-atom scattering length. The trapping 
potential experienced by a dimer has twice the value of 
the corresponding potential experienced by an atom at 
the same position. The position dependence of the num- 
ber densities and trap potentials are not shown explicitly. 
Since various recombination rates are sufliciently small 
[2,5,6,13] one can consider that the number of fermions 
and bosons are conserved sc^parately and one can account 
for the corresponding slow changes later-on when and if 
needed. As it was mentioned a few times above and it 



was shown in Ref. [5], the fermions, in the presence of 
superfluid bosons, become superfluid at sufficiently low 
temperatures. The gain in the energy density arising 
from the appearance of superfluid correlations could in 
principle be included in the above energy density. How- 
ever, this gain is very modest and is not expected to lead 
to any major changes of the spatial profile of this sys- 
tem. This gain in energy is called condensation energy 
and per particle is given by SA'^/Sep <C £f, where £ir is 
the Fermi energy. 

In the case of a large number of atoms one can find the 
atom and dimer number densities with excellent accuracy 
in the Thomas-Fermi approximation [16,17], namely form 
the equations: 
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I have absorbed in the definition of the chemical potential 
fib the energy of the dimer S2 ■ These two equations have 
been obtained by assuming that both Uf^b 7^ 0. If either 
one of the number densities is vanishing, then one has to 
solve instead one of the following two equations: 

^ +Uffnf = (If — V, if Uf > and nb = (6) 

nb = ^^^77—^, if nb > and n/ = (7) 

Ubb 

As usual in the Thomas-Fermi approximation, in the 
region of space where densities arc exponentially small 
(classically forbidden regions), the corresponding densi- 
ties are set equal to zero. In principle one can rather 
straightforwardly improve on this approximation, but 
changes in distributions are quantitatively minor. Since 
1 — 2Ufb/Ubb ~ —5, see Eq. (4), the effective trap poten- 
tial experienced by fermions changes its character in the 
presence of dimers and the fermions are pushed outward. 
In the spatial region where both n/ 7^ and nb ^ 
the fermion number density increases outwardly. By an- 
alyzing various possibilities one is quickly convinced that 
besides the rather rare cases, when there are either only 
fermions or only bosons, a finite cloud consisting of both 
Fermi atoms and Bose dimers has fermions only in the 
outside layer. 

I shall describe now the profiles of the atom and dimer 
number densities, by constructing them from point to 
point. I assume that the trap potential is increasing 
monotonically from its "center" outward. The "center" 
of the trap is defined as the point where the trapping 
potential has a minimum. The shape of the potential 
could otherwise be arbitrary. I shall fix the origin of the 
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coordinate system in the "center" of the trap and con- 
struct the number densities along an arbitrary straight 
Una. Let me denote by R3 the farthest point from origin 
along this line where the atom density vanishes, that is 
where n/(i23) = 0. At this point 

M/ = V{Rs), (8) 

and one can thus fix the value of the atom chemical po- 
tential. According to Eq. (8) the atom number density 
vanishes along an equipotential surface. Starting from 
the point /?3 one moves inward and point by point solves 
Eq. (6) and thus determine the atom number density 
along this chosen line. Now I shall assume that there are 
dimers present in the system as well, and that their num- 
ber density vanishes along the same line at the farthest 
point from origin R2 {R2 < Ra)- According to Eq. (5) 
at this point the following condition is fulfilled 

/ib = 2V{R2) + Ufi,nf{R2), (9) 

which thus determines the value of the dimer chemical 
potential. Again, the locus of points where nf,(i?2) = 
is another equipotential surface of the trapping poten- 
tial. Prom this point inward, one has now to solve the 
equations (4) and (5). One can then find a point Ri 
{Ri < R2) along this chosen line, where nf{Ri) = 0. 
Such a point might or might not exist, and its presence 
is determined by the particular values of the chemical po- 
tentials fif^b, and therefore by the particular total number 
of atoms and dimers. Naturally, if the number of dimers 
is sufficiently large, such a point where n/(ili) = 0, 
exists. If such a point Ri exists, from that point on in- 
ward the atom number density vanishes identically and 
the dimer number density is determined from Eq. (7) 
alone. Moreover, the dimer number density continues to 
increase further on from the point Ri until it reaches the 
origin. The equipotential surface, which separates the 
bosons only region at the center of the trap, from the 
region with both fermions and bosons can be determined 
from the equation: 

Ufbf^b - UbblJif = {2Ufb - Ubb)V{Ri). (10) 

One could have generated the same atom and dimer den- 
sity profiles by starting from the origin instead. In that 
case one chooses a particular value for the dimer chem- 
ical potential /ib. One has however at that point to de- 
cide whether at the origin one has atoms as well or not. 
Depending on the presence of atoms at the origin, one 
proceeds in the reverse order and generates both atom 
and dimer number densities. 

Once the number densities have been determined along 
one chosen line, one simply extends the values of the 
number densities along equipotential surfaces in the en- 
tire space. This overall onion structure is very reminis- 
cent of the structures one has inside the Earth or neutron 



stars, with a core of one nature, followed by a mantle, 
which has a structure of a different nature, and even- 
tually by a crust with a different structure as well. In 
the present case however all these three substructures 
are fluid. The extension of the present analysis to finite 
temperatures is straightforward, and one has to consider 
the free energy density, instead of the energy density, 
by adding the corresponding contributions arising from 
the fermionic and bosonic entropies as well as allowing 
for the bosons to have a non-condensed fraction. I shall 
limit the analysis here to the zero temperature case only. 
It is expected that with increasing temperature, first the 
fermionic superfluidity of the mantle is lost and subse- 
quently the superfluidity of the bosonic core. 

Some rather simple relations between these various 
radii can be obtained by taking into account the fact 
that the system is dilute, thus n/a"^ <C 1 and Ufba^ <C 1. 
One then easily obtains from Eqs. (8 - 10) that the three 
radii discussed above in a harmonic trap are related (with 
a 2% error) 

6i?i = Rl + 5Rl. (11) 

Even though terms like U jbUj have been dropped when 
this relation was obtained, the interaction is responsi- 
ble for the specific numbers appearing here, as well as 
for the existence of this relation. In particular, this re- 
lations shows that the inner core disappears as soon as 
R2 < i?3/\/6. The atom and dimer number density dis- 
tribution are then given in a spherical trap by, see Fig. 
1, 

nf = ^2{^)\^-''Rl-^-''^+^'f"^ (12) 
if Ri<r< R2, 

n/ = ^(^)'(i?3^-r2)^^ if R2<r<Rs, (13) 

The magnitude, but not the shape, of the dimer distri- 
bution depends on the particular value of the scattering 
length. And since it has a very simple shape it is not 
shown in this figure. The specific numbers here are uni- 
versal in the regime of large positive scattering lengths, 
when a » ro = {CQm/h^Y^^ . Here w is the trap fre- 
quency and Cg is the strength of the atom-atom van der 
Waals interaction. The atom distribution is thus essen- 
tially independent of the particular value of the scattering 
length a, as long as this is large and the system is dilute. 
One can easily derive analytical formulas for the num- 
ber of atoms and dimers as functions of the radii -Ri,2,3, 
even for anisotropic traps. Once the entire system is in 
the quantum degenerate regime the temperature depen- 
dence of the number density profiles are rather weak [16]. 
The ratio of dimers to atoms scales as a function of w, a 
and cloud radii as follows 
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FIG. 1. The atom density profile of an atom-dimer mixture 
in the dilute limit, when R2 = 0.81 R3 and where the number 
density was normalized to one. In this case about 16% of the 
atoms are in the "mantle." 

-^7 OC -f.. (15) 

Since the scattering length a is the smahest length scale 
typically the number of dimers is larger than the number 
of atoms, but this is not true always. 

Besides the rather intricate structure of this new quan- 
tum object, one can naturally expect that the spectrum 
of collective states is going to have a nontrivial charac- 
ter, arguably different from the spectrum of BEC sys- 
tems studied so far [16,18]. As noted above, the core 
and the mantle of this new quantum object are super- 
fluid, the first as a Bose system and the latter as mixed 
Bose-Fermi superfluids. As the most notable feature of a 
superfluid system is the possibility to excite a superflow, 
it would be extremely interesting to study the charac- 
ter of one or more vortex states. Unlike Bose [18,19] or 
Fermi [20] clouds, since the crust here is a normal Fermi 
fluid, a vortex will have to be limited to either the core 
and mantle or the mantle alone, and not extend into the 
Fermi normal crust (apart from the proximity effects due 
to the Andreev reflection mentioned above). It is likely 
that the Bose and Fermi superfluids can sustain largely 
independent superflows and the presence of a vortex in 
one component would not necessarily lead to the appear- 
ance of a vortex in the other component as well. The 
rotational properties of this object are clearly nontrivial, 
as one can expect both rotational and irrotational flows 
at the same time. The kinetics and thermalization of an 
object with such a complex structure is going to provide 
a interesting challenge to unravel, both experimentally 
and theoretically. However, unlike stars and planets, we 
should be able to study these processes under a variety of 
conditions and be able to finely tune various properties 
of these "stars in a bottle." 

In summary, it was shown that an atomic cloud consist- 
ing of fermionic atoms in a trap, at sufficiently low tem- 
peratures and when the atom-atom scattering length is 



positive and large, will have an onion-like structure, with 
a superfluid bosonic core made of atomic dimers, a man- 
tle consisting of a mixture of superfluid bosonic dimers 
and fermionic atoms, all surrounded by an outer layer 
consisting of fermions alone, forming a normal Fermi gas. 
The existence of this particular structure is due to a large 
extent to the fact that various couplings "conspired" to 
acquire the specific values quoted here, see Eqs. (3-5). 
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